Full Counting Statistics in Strongly Interacting Systems: Non-Markovian Effects 
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We present a theory of full counting statistics for electron transport through interacting electron 
systems with non-Markovian dynamics. We illustrate our approach for transport through a single- 
level quantum dot and a metallic single-electron transistor to second order in the tunnel-coupling 
strength, and discuss under which circumstances non-Markovian effects appear in the transport 
properties. 
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The study of current fluctuations in mesoscopic sys- 
tems has become an intense field of research, since it 
allows to access information about electron correlations 
that is not contained in the average current. The phe- 
nomenon of shot noise 0, S 13 > that dominates the cur- 
rent noise at low temperatures, has been investigated the- 
oretically and experimentally in various contexts. Fur- 
ther understanding in electron transport can be gained 
from the study of higher moments Q which can be conve- 
niently extracted from the study of full counting statistics 
(FCS) Up to date FCS has been studied in a vari- 

ety of cases. Exaroples include normal-superconductor 
hybrid structures j3|, superconducting weak links 0, 
tunnel junctions 0, chaotic cavities [iffl. e ntangled elec- 
trons [ij and spin-correlated systems [iJI- Schemes for 
an experimental measurement of FCS has been put for- 
ward in Refs. 

Electron-electron interaction may strongly influence 
quantum transport. The effect of electron correlation 
on FCS has been considered so far in the case of a 
weakly interactin g rn esoscopic conductor |l4| . almost 
open dot s [liL ImIt^ , charge pumping 0] and charge 
shuttles [1^. The FCS for Coulomb-blockade devices 
has been analyzed 0, IIS 1^ in the framework of a 
Born-Markov master equation approach. The aim of the 
present paper is to extend this idea to a obtain a general 
theory of full counting statistics for strongly interacting 
systems with non-Markovian behavior. In particular, we 
formulate a perturbative non-Markovian expansion that 
allows for a systematic study of the relative importance of 
non-Markovian corrections. We demonstrate for the ex- 
ample of a single-level quantum dot with strong Coulomb 
interaction that non-Markovian effects become increas- 
ingly important for higher moments of the current fluc- 
tuations. 

Full counting statistics. - Full information about all 
transport properties of a given system is contained in 
the probability distribution P{N,t) that N charges have 
passed through the system after time t. This distribution 
function is related to the cumulant generating function 



Six) by 



5(X) = -In 



J2 e'^^P{N,t) 



(1) 



where x is the counting field. All moments of the current 
can be obtained from the cumulant generating function 
by performing derivatives with respect to the counting 
field ((/))„ = -(-i)"(e"/i)5;^S'(x)lx=o- The first four 
moments are the average current, the (zero- frequency) 
current noise, the skewness, and the kurtosis. 

In this work we consider systems with strong local in- 
teractions, such as electrons in a quantum dot, that are 
coupled to a reservoir of noninteracting degrees of free- 
dom. In these situations transport properties can often 
be described in terms of few (local) degrees of freedom 
(the charge of the quantum dot in the previous example). 
It is then convenient to integrate out the noninteracting 
degrees of freedom to arrive at an effective description of 
the reduced system only. Let p'" be the vector of prob- 
abilities to find the system in the corresponding state at 
the initial time t = Q. The time evolution of the system 
is described by a generalized master equation. 



dt 



p(7V,t) 



E 

Ar'=-c 



dt' W{N - N',t,t') -piN^t') , 



(2) 

where p{N,t) is the vector of dot occupation probabil- 
ities under the condition that N electrons have passed 
the system. The cumulant generating function is given 
by Eq. ^ where P{N,t) ^ e'^ • p{N,t) with = 
(1,1,...,!). The matrix W(A^ - N', t, t') describes tran- 
sitions during which N — N' electrons are transferred. 
The counting field x is introduced by Fourier transform- 
ing the master equation p(x7*) = 6xp(«-/Vx)p(-^7 ^) 
and W(x, t, t') = Ejv eM^Nx)W{N, t, t'). 

In general, the kernels W(x,t,i') are nonlocal in time 
as a consequence of having integrated out the reservoir 
degrees of freedom. We consider the case in which there 
is no explicit time dependence of the systems parame- 
ters, so that W(x, t — t') can be Laplace transformed. 



2 



W{x,z) = dt exp{-zt)W{x,t). Usually, the dy- 
namics of the system is characterized by a time scale of 
the individual transitions. We assume that W(x,t) de- 
cays faster than any power of t, such that the derivatives 
9"W(x, exist for all n. 

In the Markovian limit, W(x,t) ~ S{t), the system 
at time t depends only on the state of the system at 
the same time t. The master equation Eq. (O is, then, 
solved by p(x, t) = exp[W(x, z = 0)t] ■ p'", i.e., only the 
z — component of the transitions W(x, z) are taken 
into account. 

Our goal is to describe stationary transport properties 
in the presence of a memory of the system, described by 
the full z-dependence of W(x, z). To solve the master 
equation without making use of the Markovian approxi- 
mation, we switch to the Laplace representation. 



zpix, z) - p'" = W(x, z) ■ p(x, z) , 



(3) 



which is solved by p(x, -z) = ^^=q[^{Xi z)]^ / z^^^p™ . 
By assuming that the kernel W(x, t) decays in time faster 
than any power law, we can define the Taylor series 
[Wix,z)r = E™=o5r[W(x,^)]"L^o2"/m! and sub- 
stitute it in the previous solution for p(x, z). The long- 
time behavior of p(Xj z) is determined by its poles in z. 
Performing the inverse Laplace transformation we get 



P(x,i)-E 



([W(x, 



]n ■W{x,z)t\ 



n=0 



z=0+ 



•P 



(4) 



for large t. To proceed, we perform a spectral decom- 
position of the matrix W(x, z). For physical reasonable 
systems, all eigenvalues have a negative real part. As 
a consequence of the exponential function in Eq. Q), 
the long-time behavior will be dominated by the eigen- 
value A(x, z) with the smallest absolute value of the 
real part. Let qo and po be the corresponding left 
and right eigenvectors, Qq • W(x, — \{x, z)qQ , and 
W(x, z) ■ Po = A(x, 2)po- Unitarity in the absence of 
counting fields implies A(0, z) — for all z. 

The cumulant generating function becomes 
Six) = - In [En=o (A"e^*+^)] , with 
fJ,{x,z) = ln[(e^ • Po)(qo ■ Pin)]- By perform- 
ing the time derivative and making use of the 
relation ^^=0 («^'"+') / ^^=0 («^") 
En=o /Er=o that holds for arbi- 

trary functions a and b, we arrive at the central result of 
this Letter, 



Six) 



j:Zoi!dn>^Hx,z)] 



t. 
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The cumulant generating function depends only on 
the eigenvalue A(x, 2). This result can be used as a 
starting point for a non-Markovian expansion. Six) ~ 
J2'^=o Snix) where S'„(x) contains n derivatives with re- 
spect to z applied to n -f- 1 factors of A. While S'o(x) 



describes the Markovian limit, Snix) is the n-th non- 
Markovian correction. 

Perturbative non-Markovian expansion. - For many 
systems there is a small parameter which allows for a 
perturbative analysis of all transport properties. In the 
examples to be discussed below this will be the tunnel- 
coupling strength between the leads and the interact- 
ing region (quantum dot or metallic island). Then, 
A(x, z) — X^i^i '^'■'•'(X; ^)! where the superscript (i) indi- 
cates the order in the small parameter. The lowest-order 
transport properties are derived from the lowest-order cu- 
mulant generatin g fu nction 5'^^^(x) = — iA*^^) (x, 2;)|z=o+ : 
as found in Ref. [19|. Non-Markovian corrections, sig- 
naled by derivatives 5*A(x, z), do not enter in this limit. 
The highest derivative that enters in the evaluation of 
the n-th moment in m-th order perturbation theory is 
9^ ^(x, 2;)|z=o with k = min{n, m} — 1. As a conse- 
quence, non-Markovian behavior is probed only in the 
second or higher moment combined with second or higher 
order in perturbation theory. The second-order con- 
tribution for example reads S'^'^\x) = ^A^^'^^ix^z) + 
A'^^)(x, 2;)9zA^^^(x, 2)]z=o+- The appearance of these 
derivatives in the noise of second-order transport through 
quantum dots has been also found in Ref. |23. In the re- 
maining part of this Letter, we illustrate our approach 
with two examples. We calculate the cumulant generat- 
ing function for second-order transport through a single- 
level quantum dot and through a metallic single-electron 
transistor in the presence of strong Coulomb interaction. 

Single level QD.- The single-level quantum dot is de- 
scribed by the Hamiltonian, H = H]^ + + Hi^ + H^. 
The electrons in the noninteracting left and right leads 
are represented by Hi^ and iJ^, respectively, Hd = 
e^^cj^Co- -l- Uuiui describes the dot with level energy 
e and charging energy U for double occupation. Tun- 
neling is modeled by H^.r = J^a^rO-lka'^'^ + with 
r = R, L, where we assume the tunnel matrix element 
tr to be independent of momentum k and spin a. The 
tunnel-coupling strength is characterized by the intrinsic 
linewidth F = Fl -l- Fr with F^ = 27rpr|^rP where pr is 
the density of states in the leads. An asymmetry of the 
tunnel couplings is parametrised by 7 = 4FlFr/F^. 

To derive the kernels W of the generalized master 
equation, we make use of a diagrammatic real-time tech- 
nique [23| for the time evolution of the reduced density 
matrix formulated on a Keldysh contour. We introduce 
counting fields Xr for tunneling through barrier r into 
lead r by the replacement tr — > exp(iXr) for tunnel 
vertices on the upper and tr — > tr exp(— ixr) on the lower 
branch of the Keldysh contour with xl = ~XR — x/2- 

We consider the limit U — *■ 00, in which double occu- 
pancy of the dot is prohibited [2J| , and obtain 



s^'\x)- 



trfie) 
2h 



1 - 



27Efc/wW(e''=^-l) 



[fi^)Y 



(6) 
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with = [1 - /lH]/r(^), /(_)(c^) = /l(u;)[1 - 

/r(w)], and f{uj) = Er^rll + /r(w)]/r, where /^(w) is 
the Fermi function for lead r. This result was previ- 
ously obtained in |l9l |. The second-order contribution, 

5'(2)(x) = 5'^^t(x) + S^elix), consists of two terms. The 
first one, 



0.20 



(e^'^^-l) ldu;f^k)i^)R{u;-e), (7) 



k=± 



with Ricj) = Re[l/(a;-|-iO'*')^], describes cotunneling pro- 
cesses 25j, and is in agreement with previous work about 
noise in cotunneling regime [2^ . The counting- field de- 
pendence corresponds to a bidirectional Poisson statistics 
of a single barrier where the transition rates are substi- 
tuted by the cotunneling rates of the quantum dot. How- 
ever, S^'^\x) contains a second contribution. 



s^^lix) = dJs<^'\x)M^ie)] 



(8) 



with a{e) = - X;^(r^/27r) / dLoU{uj)/[uj - e + iO+). In 
the previous formulas an high energy cut-off E^, of the 
order of charging energy, has to be introduced in order 
to cure spurious divergences related to the fact that we 
restricted the charge states to 0,1 ^3- The contribu- 
tion S'rcn(x) is also of second order in the tunnel-coupling 
strength but obeys the same statistics as the first-order 
(sequential-tunneling) result, Eq. ©. This suggests that 
there are two different types of second-order contribu- 
tions to transport. In addition to the usual cotunneling 
processes, there are corrections to sequential tunneling 
due to quantum-fluctuation induced renormalization of 
the system parameters. From the form of Eq. © we 
deduce a renormalization of level position and coupling 
strength given by e = e-l-Recr(e) and F = r[l-|-i9eRe(7(e)]. 

In Fig.^we plot the flrst four moments (current, noise, 
skewness, and kurtosis) as a function of level position 
e. The solid lines represent the full first- plus second- 
order result, as compared to the first-order contribution 
(dashed fine). 

The relative importance of the non-Markovian con- 
tributions is illustrated in Fig. [3 While for the cur- 
rent only Markovian contributions enter (see discussion 
above), non-Markovian corrections become increasingly 
important for higher moments. 

Metallic QD- A similar analysis can be performed 
for a metallic single-electron transistor, which accom- 
modates a continuum of states on the dot and in- 
cludes a large number of transverse channels. Follow- 
ing the notation of Refs. I27I we characterize the tunnel- 
coupling strength by the dimensionless conductance — 
h/ (Aire'^Rr) where Rr is the resistance of barrier r = L, R. 
We concentrate on the low-temperature regime, in which 
only two charge states of the metallic island have to be 
taken into account. This requires, again, the introduc- 
tion of an high-energy cut-off to regularize the integrals. 
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FIG. 1: Current (a), noise (b), skewness (c), and kurtosis 
(d) for a single-level quantum dot with (/ — > 00 as a function 
of the level position e/eVsd for first (dashed lines) and first- 
plus second-order (solid lines) in tunneling with eVsd = (mr — 
/^l). Other parameters are 7 = 1, F = 3 ■ 10~^e\4d, ksT = 
10~^eV^d and the high-energy cutoff Ec = 10eX4d. 



The difference of the charging energy between them is 
denoted by A. We obtain for the first order of cumulant 
generating function 



s^'Hx) 



tTra{A) 



\ 



a(fc)(A)(e»fcx 



1) 



fc=± 



a(A) 



(9) 

where we have used the definitions a{u!) — J2r=L r o^ri^) 
and Q!(±)(a;) — a^{u!)a^{uj) with a^{uj) = ar(uj)f^{uj) 
where {uj) is the Fermi functions of lead r and ar{uj) = 
Qq coth[/3((jj — /ir)/2]. The second-order term is given by 

S^'^Kx) = 'S'cot(x) + Srcl{x)i as the case of the single- 
level quantum dot. Again, there is a contribution due to 
cotunneling. 



siolix) 



and a term 



2TTt 



fc=± 



ikx 



1) /dwa(fc)(w)i?(w - A), 
(10) 



si'2ix)^dA\s^'Hx)M^m 



(11) 



associated with sequential-tunneling processes with 
renormalized system parameters, where cr(A) = 
— J2rl rf'^ar ('^)/(w— A-|-iO+). This interpretation of the 
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FIG. 2: Relative contribution of the non-Markovian part 
kki))^ ~ to the n-th moment {{/))„ in first- plus 

second-order in tunneling. The full, dashed, and dot-dashed 
line corresponds to noise in — 2), skewness (n = 3), and 
kurtosis (n — 4) for the same parameters as in Fig. 



different types of second-order contributions is consistent 
with the analysis of the second-order current 28] and of 
the FCS within a drone Majorana fermion representation 

Conclusions. ~ We present a theory of FCS for inter- 
acting systems with non-Markovian dynamics. A general 
expression for the cumulant generating function is de- 
rived that provides the starting point for a perturbative 
non-Markovian expansion. As examples we study trans- 
port through a single-level quantum dot and a metallic 
single-electron transistor to second order in the tunnel- 
coupling strength. From our formulation we could iden- 
tify two different types of contributions to second-order 
transport, namely cotunneling and corrections to sequen- 
tial tunneling due to renormalization of the system pa- 
rameters. Furthermore, we demonstrate the increasing 
importance of non-Markovian effects for higher moments 
and higher orders in the tunnel- coupling strength. 
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